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EDITORS' PREFACE. 



The volume called Higher Mathematics, the first edition 
of which was published in 1896, contained eleven chapters by 
eleven authors, each chapter being independent of the others, 
but all supposing the reader to have at least a mathematical 
training equivalent to that given in classical and engineering 
colleges. The publication of that volume is now discontinued 
and the chapters are issued in separate form. In these reissues 
it will generally be found that the monographs are enlarged 
by additional articles or appendices which either amplify the 
former presentation or record recent advances. This plan of 
publication has been arranged in order to meet the demand of 
teachers and the convenience of classes, but it is also thought 
that it may prove advantageous to readers in special lines of 
mathematical literature. 

It is the intention of the publishers and editors to add other 
monographs to the series from time to time, if the call for the 
same seems to warrant it. Among the topics which are imder 
consideration are those of elliptic functions, the theory of num- 
bers, the group theory, the calculus of variations, and non- 
Euclidean geometry; possibly also monographs on branches of 
astronomy, mechanics, and mathematical physics may be included. 
It is the hope of the editors that this form of publication may 
tend to promote mathematical study and research over a wider 
field than that which the former volume has occupied. 

December, 1905. 
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AUTHOR'S PREFACE. 



In republishing this short treatise in book form the author 
solicits criticism but offers no apology. The type of the book 
he has sought to imitate is that shown in the " mathematical 
tracts " of the late Sir George B. Airy. The brevity and the con- 
crete illustrations of these " tracts " have served very effectively 
in introducing students to a number of the more difficult fields 
of applied mathematics; and it is hoped that this treatise will 
serve a similar end. 

The theory of probability and the theory of errors now con- 
stitute a formidable body of knowledge of great mathematical 
interest and of great practical importance. Though developed 
largely through applications to the more precise sciences of as- 
tronomy, geodesy, and physics, their range of applicability extends 
to all of the sciences; and they are plainly destined to play an 
increasingly important r6le in the development and in the appli- 
cations of the sciences of the future. Hence their study is not 
only a commendable element in a liberal education, but some 
knowledge of them is essential to a correct understanding of 
daily events. 

No special novelty of presentation is claimed for this work; 
but the reader may find it advantageous to know that a definite 
plan has been followed. This plan consists in presenting each 
principle, first, by means of a simple, concrete example; passing, 
secondly, to a general statement by means of a formula; and, 
thirdly, illustrating applications of the formula by concrete 
examples. Great pains have been taken also to secure clear and 
correct statements of fundamental facts. If these latter are 
duly understood, the student needs little additional aid; if 
they are not duly understood, no amount of aid will forward him. 

The passage from the elementary concrete to the advanced 
abstract may appear to be abrupt to the reader in some cases. 
It is hoped, however, that any large gaps may be easily bridged 
and that any serious difficulties may be easily overcome by means 
of the references given to the literature of the subject. In any 
event the student will find that in this, as in all of the more ardu- 
ous sciences, his greatest pleasure and his highest discipline will 
come from bridging such gaps and from surmounting such diffi- 
culties. 

Washington, D. C, December, 1905. 
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PROBABILITY AND THEORY OF ERRORS. 



ART. 1. Introduction. 

It is a curious circumstance that a science so profoundly 
mathematical as the theory of probability should have origi- 
nated in the games of chance which occupy the thoughtless and 
the profligate.* That such is the case is sufficiently attested 
by the fact that much of the terminology of the science and 
many of its familiar illustrations are drawn directly from the 
vocabulary and the paraphernalia of the gambler and the trick- 
ster. It is somewhat surprising, also, considering the antiquity 
of games of chance, that formal reasoning on the simpler 
questions in probability did not begin before the time of Pascal 
and Fermat. Pascal was led to consider the subject during the 
year 1654 through a problem proposed to him by the Chevalier 
de M^r^, a reputed gambler.f The problem in question is 
known as the problem of points and may be stated as follows : 
two players need each a given number of points to win at a 
certain stage of their game ; if they stop at this stage, how should 
the stakes be divided? Pascal corresponded with his friend 
Fermat on this question ; and it appears that the letters which 
passed between them contained the earliest distinct formulation 
of principles falling within the theory of probability. These 

* The historical facts referred to in this article are drawn mostly from Tod- 
hunter's History of the Mathematical Theory of Probability from the time of 
Pascal to that of Laplace (Cambridge and London, 1865). 

f ** Un probl^me relatif aux jeux de hasard, propos6 iun austere jans6niste 
par un homme du monde, a 6t6 Torigine du calcul des probabilit6s. " Poisson, 
Recherches sur la Probability des Jugements (Paris, 1837). 
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acute thinkers, however, accomplished little more than a correct 
start in the science. Each seemed to rest content at the time 
with the approbation of the other. Pascal soon renounced 
such mundane studies altogether ; Fermat had only the scant 
leisure of a life busy with affairs to devote to mathematics; 
and both died soon after the epoch in question, — Pascal in 
1662, and Fermat in 1665. 

A subject which had attracted the attention of such dis- 
tinguished mathematicians could not fail to excite the interest 
of their contemporaries and successors. Amongst the former 
Huygens is the most noted. He has the honor of publishing 
.the first treatise* on the subject. It contains only fourteen 
propositions and is devoted entirely to games of chance, but it 
gave the best account of the theory down to the beginning of 
the eighteenth century, when it was superseded by the more elab- 
orate works of James Bernoulli,t Montmort,:]: and De Moivre.§ 
Through the labors of the latter authors the mathematical 
theory of probability was greatly extended. They attacked, 
quite successfully in the main, the most difficult problems ; 
and great credit is due them for the energy and ability dis- 
played in developing a science which seemed at the time to 
have no higher aim than intellectual diversion.! Their names, 
undoubtedly, with one exception, that of Laplace, are the most 
important in the history of probability. 

Since the beginning of the eighteenth century almost every 
mathematician of note has been a contributor to or an expos- 
itor of the theory of probability. Nicolas, Daniel, and John 
Bernoulli, Simpson, Euler, d'Alembert, Bayes, Lagrange, Lam- 
bert, Condorcet, and Laplace are the principal names which 
figure in the history of the subject during the hundred years 

*De Ratiociniis in Ludo Ales, 1657. 

t Ars Conjectandi, 1713. 

X Essai d'Analyse sur les Jeux de Hazards, 1708. 

§The Doctrine of Chances, 1718. 

|Todhunter says of Montmort, for example, ** In 1708 he published his 
work on Chances, where with the courage of Columbus he revealed anew world 
to Mathematicians." 
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ending with the first quarter of the nineteenth century. Of 
contributions from this briUiant array of mathematical talent, 
the Th^orie Analytique des Probabilitts of Laplace is by far 
the most profound and comprehensive. It is, like his M&- 
canique Celeste in dynamical astronomy, still the most elabo- 
rate treatise on the subject. An idea of the grand scale of the 
work in its present form* may be gained by the facts that the 
non-mathematical introduction! covers about one hundred and 
fifty quarto pages; and that, in spite of the extraordinary 
brevity of mathematical language, the pure theory and its ac- 
cessories and applications require about six hundred and fifty 
pages. 

From the epoch of Laplace down to the present time the 
extensions of the science have been most noteworthy in the 
fields of practical applications, as in the adjustment of obser- 
vations, and in problems of ijisurance, statistics, etc. Amongst 
the most important of the pioneers in these fields should 
be mentioned Poisson, Gauss, Bessel, and De Morgan. Nu- 
merous authors, also, have done much to simplify one or an- 
other branch of the subject and tlms bring it within the range 
of elementary presentation. The fundamental principles of 
the theory are, indeed, now accessible in the best text-books 
on algebra : and there are many excellent treatises on the pure 
theory and its various applications. 

Of all the applications of the doctrine of probability none 
is of greater utility than the theory of errors. In astronomy, 
geodesy, physics, and chemistry, as in every science which at- 
tains precision in measuring, weighing, and computing, a 
knowledge of the theory of errors is indispensable. By the aid 
of this theory the exact sciences have made great progress dur- 
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ing the nineteenth century, not only in the actual determination 
of the constants of nature, but also in the fixation of clear 
ideas as to the possibilities of future conquests in the same di- 
rection. Nothing, for example, is more satisfactory and in- 
structive in the history of science than the success with which 
the unique method of least squares has been applied to the 
problems presented by the earth and the other members of the 
solar system. So great, in fact, are the practical value and 
theoretical importance of the method of least squares, that it is 
frequently mistaken for the whole theory of errors, and is 
sometimes regarded as embodying the major part of the doc- 
trine of probability itself. 

As may be inferred from this brief sketch, the theory of 
probability and its more important applications now constitute 
an extensive body of mathematical principles and precepts. 
Obviously, therefore, it will be impossible within the limits of 
a single condensed monograph to do more than give an out- 
line of the salient features of the subject. It is hoped, how- 
ever, in accordance with the general plan of the volume, that 
such outline will prove suggestive and helpful to those who 
may come to the science for the first time, and also to those 
who, while somewhat familiar with the difficulties to be over- 
come, have not acquired a working knowledge of the subject. 
Effort has been made especially to clear up the difficulties of 
the theory of errors by presenting a somewhat broader view of 
the elements of the subject than is found in the standard 
treatises, which confine attention almost exclusively to the 
method of least squares. This chapter stops short of that 
method, and seeks to supply those phases of the theory which 
are either notably lacking or notably erroneous in works 
hitherto published. It is believed, also, that the elements here 
presented are essential to an adequate understanding of the 
well-worked domain of least squares.* 

*The author has given a brief but comprehensive statement of the method 
of least squares in the volume of Geographical Tables published by the Smith- 
sonion Institution, 1894. 
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Art. 2. Permutations. 
The formulas and results of the theory of permutations 
and combinations are often needed for the statement and so- 
lution of problems in probabilities. This theory is now to be 
found in most works on algebra, and it will therefore suffice 
here to state the principal formulas and illustrate their mean- 
ing by a few numerical examples. 

The number of permutations of « things taken r in a group 
is expressed by the formula 

(»), = „(« -lX»- 2)... (»-r+ I). (1) 

Thus, to illustrate, the number of ways the four letters a, b, 
c, d can be arranged in groups of two is 4 . 3 = 1 2, and the groups 
are 

dc. 



ah, ba, ac, ca, ad, da, be, cb, bd 


db, ed. 


Similarly, the formula gives for 




n = 3 aiidr = 2, (3),= 3-2 


= 6, 


«= 7 " '■ = 3. (7). = 7-6-S 


= 210 



w = lo " r = 6, (io),=: 10.9.8.7.6, 5 = 151200. 
The results which follow from equation (i) when n and r 
do not exceed lo each are embodied in the following table: 
. Values of Permui 
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CS64100 


986409 


,0.600 


13699 


1056 


325 


64 


'5 



The use of this table is obvious. Thus, the number of per- 
mutations of eight things in groups of five each is found in the 
fifth line of the column headed with the number 8. It will be 
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noticed that the last two numbers in each column (excepting 
that headed with i) are the same. This accords with the for- 
mula, which gives for the number of permutations of n things 
in groups of n the same value as for n things in groups of 
{n — i). It will also be remarked that the last number in each 
column of the table is the factorial, «!, of the number n at the 
head of the column. For example, in the column under 7, the 
last number is 5040 =1.2.3.4.5.6.7 = 7!. 

The total number of permutations of n things taken singly, 
in groups of two, three, etc., is found by summing the numbers 
given by equation (i) for all values of r from I to n. Calling 
this total or sum Spy it will be given by 

Sp=:S{n)^ (2) 

To illustrate, suppose « = 3, and, to fix the ideas, let the 
three things be the three digits I, 2, 3. Then from the above 

table it is seen that 5^ = 3 + ^ + ^ = '5 5 ^^* ^^^^ ^^^ number 
of numbers (all different) which can be formed from those dig- 
its is fifteen. These numbers are i, 2, 3 ; 12, 13, 21, 23, 31, 32 ; 
123, 132, 213, 231, 312, 321. 

The values of 5^ for « = i, 2, ... 10 are given under the 
corresponding columns of the above table. But when n is 
large the direct summation indicated by (2) is tedious, if not 
impracticable. Hence a more convenient formula is desirable. 
To get this, observe that (i) may be written 

if r IS restricted to integer values between i and {n — i), both 
inclusive. This suffices to give all terms which appear in the 
right-hand member of (2), since the number of permutations 
for r = (« ~ i) is the same as for r = «. Hence it appears 
that 

^ ' I ' I . 2 ' (« — i)! 



(«~ I) 

= « !( I H h h . . . 7 Ti)' 

\ ^ I ' 1 . 2 ' (» — i) !/ 
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But as n increases, the series by which n !. is here multiplied 
approximates rapidly towards the base of natural logarithms; 
that is, towards 

e = 2.7 1 828 1 8 +» log ^ = 04342945. 

Hence for large values of n 

Sp= nle, approximately.* (3) 

To get an idea of the degree of approximation of (3), sup- 
pose n = g. Then the computation runs thus (see values in 
the above table) : 

log 

9! = 362880 5-5597630 
e 0.4342945 

9! ^ = 986410 5-9940575 
Sp = 986409 by equation (2). 

The error in this case is thus seen to be only one unit, or 
about one-millionth of Sp,\ 

Prob. I. Tabulate a list of the numbers of three figures each 
which can be formed from the first five digits i, . . . 5. How many 
numbers can be formed from the nine digits ? 

Prob. 2. Is Sp always an odd number for n odd ? Observe 
values of •S^ in the table above. 

Art. 3. Combinations. 

In permutations attention is given to the order of arrange- 
ment of the things considered. In combinations no regard is 
paid to the order of arrangement. Thus, the permutations of 
the letters ^, by c, d in groups of three are 

(abc) {abd) bac bad acb {acd) cab cad 
adb adc dab dac bca {bed) cba cbd 
bda bdc dba dbc cda cdb dca deb 

* See Art. 6 for a formula for computing n\ when /« is a large number. 

f When large numbers are to be dealt with, equations (i)' and (3) are easily 
managed by logarithms, especially if a table of values of log («!) is available. 
Such tables are given to six places in De Morgan's treatise on Probability in 
the Encyclopaedia Metropolitana, and to five places in Shortrede's Tables 
(Vol, I, 1849). 
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But if the order of arrangement is ignored all of these are 
seen to be repetitions of the groups enclosed in parentheses, 
namely, {abc), {ctbd), (acd)^ {bed). Hence in this case out of 
twenty-four permutations there are only four combinations. 

A general formula for computing the number of combina- 
tions of n things taken in groups of r things is easily derived. 
For the number of permutations of n things in groups of r is 
by (i) of Art. 2 

{n)r = n{ti — i){n — 2) . , . (« — r -f- i) ; 

and since each group of r things gives 1.2.3. ..r = r! P^^" 
mutations, the number of combinations must be the quotient 
of {n\ by r!. Denote this number by Cin)^. Then the gen- 
eral formula is 

, _ n{n-- !)(;/- 2).. .(;/-r+ i) 
^Wr = -\ (I) 

This formula gives, for example, in the case of the four let- 
ters tf, b, c, d taken in groups of three, as considered above, 

Multiply both numerator and denominator-of the right-hand 
member of (i) by (« — r) ! The result is 

C{n)r =>!(«!;.)!' (0' 

which shows that the number of combinations of n things in 
groups of r is the same as the number of combinations of n 
things in groups of {n — r). Thus, the number of combina- 
tions of the first ten letters a, b, c . . .j in groups of three or 

seven is 

10! 
—, — j = 120. 

3!?! 
The following table gives the values C{n)r for all values of 
n and r from I to 10. 

The mode of using this table is evident. For example, the 
number of combinations of eight things in sets of five each is 
found on the fifth line of the column headed 8 to be 56. 



combinations. 
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It will be observed that the numbers in any column show 
a maximum value when h is even and two equal maximum 
values when w Is odd. That this should be so is easily seen 
from (1)', which shows that €{»), will be a maximum for any 
value of It when r ! {« — r) ! is a minimum. For w even this is 
a minimum for r = ^n; while for « odd it has equal minimum 
values for r = J(« — i) and r = ^(« + i). Thus, 



1 of an\ = 



^ for « even, 



(2) 



The total number of combinations of « things taken singly, 
in groups of two, three, etc., is found by summing the numbers 
given by (i) for all values of r from i to « both inclusive. 
Calling thi3 total or sum 5^, 

5, = SCin\. 
The same sum will also come from (1)' by giving to r all values 
from I to (« — i), both inclusive, summing the results, and in- 
creasing their aggregate by unity. Thus by either process 



5, = « + - 



-h « + I. 
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The second member of this equation is evidently equal to 

(i + 0* "■ ^- Hence 

5, = 2C{n)r = 2- - I. (3) 

The values of S^ for values of n and r from i to lo are given 
under the corresponding columns of the above table. 

Prob. 3. How many different squads of ten men each can be 
formed from a company of 100 men ? 

Prob. 4 How many triangles are formed by six straight lines 
each of which intersects the other five ? 

Prob. 5. Examine this statement : " In dealing a pack of cards 
the number of hands, of thirteen cards each, which can be produced 
is 635 013 559 600. But in whist four hands are simultaneously held, 
and the number of distinct deals . . . would require twenty-eight 
figures to express it." * 

Prob. 6. Assuming combination always possible, and disregarding 
the question of proportions, find how many different substances 
could be produced by combining the seventy-three chemical ele- 
ments. 

Art. 4. Direct Probabilities. 

If it is known that one of two events must occur in any 
trial or instance, and that the first can occur in a ways and the 
second in d ways, all of which ways are equally likely to hap- 
pen, then the probability that the first will happen is expressed 
mathematically by the fraction a/{a-\-b)y while the probability 
that the second will happen is b/(a -|- b\ Such events are said 
to be mutually exclusive. Denote their probabilities by p and 
q respectively. Then there result 

the last equation following from the first two and being the 
mathematical expression for the certainty that one of the two 
events must happen. 

Thus, to illustrate, in tossing a coin it must give " head *' or 
" tail" ; ^? = ^ = I, and p z= q =z 1/2. Again, if an urn contain 
^ = 5 white and ^ = 8 black balls, the probability of drawing 

* Jevons, Principles of Science, New York, 1874, p. 217. 
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a white ball in one trial is / = 5/13 and that of drawing a 
black one q = 8/13. 

Similarly, if there are several mutually exclusive events 
which can occur in a, b, c , , , ways respectively, their probabil- 
ities p, q, r . . , are given by 

a b __ ^ 

(2) 
P + 9 + ^ + " •= I- 

For example, if an urn contain a=4 white, b= $ black, 
and c = 6 red balls, the probabilities of drawing a white, black, 
and red ball at a single trial are p = 4/i$, ^=5/15, and 
r = 6/15, respectively. 

Formulas (i) and (2) may be applied to a wide variety o.f 
cases, but it must suffice here to give only a few such. As a 
first illustration, consider the probability of drawing at random 
a number of three figures from the entire Hst of numbers which 
can be formed from the first seven digits. A glance at the 
table of Art. i shows that the symbols of formula (i) have in 
this case the values a= 210, and a-^-b^ 13699. Hence 
b = 13489, and / = 210/13699 ; that is, the probabiHty in ques- 
tion is about 1/65. 

Secondly, what is the probability of holding in a hand of 
whist all the cards of one suit ? Formula (i) of Art. 3 shows 
that the number of different hands of thirteen cards each which 
may be formed from a pack of fifty-two cards is 

52.51. 50... 40 ^ ^ 

..2.3. ..13 =^35013559600, 

and the probability required is the reciprocal of this number. 
The probability against this event is, therefore, very nearly 
unity. 

Thirdly, consider the probabilities presented by the case of 
an urn containing 4 white, 5 black, and 6 red balls, from which 
at a single trial three balls are to be drawn. Evidently the 
triad of balls drawn may be all white, all black, all red, partly 
white and black, partly white and red, partly black and red, or 
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one each of the white, black, and red. There are thus seven 
different probabiHties to be taken into account. The theory 
of combinations shows (see equation (i), Art. 3) that the total 
number of 

White triads = — '-^—^ — = 4 = a 

1.2.3 

Black triads = —1— =10 = ^ 

6 

^ , . , 6.5.4 

Red triads = — z — = 20 = c 



I . 


2 . 


3 


5' 


.4. 
6 


3 


6. 


■5- 


4 




6 




9 


.8. 


7 




b 




10 


■9 


.8 




6 




II 


.10.9 



White and black triads = — ^ — ( 4+10)= 70 = d 

10. 9 . 8 
White and red triads = — ^ — — ( 4-|-2o)= g6 = e 

Black and red triads — — 1— _ !-^ — (io-|-2o) = i35 =/ 

White, black, and red triads = 4.5.6 =ri20=-^ 

Sum =455 

The total number of these triads is 455, and is, as it should 
be, the number of combinations in groups of three each of the 
whole number of balls. Hence formulas (2) give the seven 
different probabilities which follow, using the initial letters 
w, by r to indicate the colors represented in a triad : 
For a triad www p = 4/455, 

did q = 10/455, 

rrr r = 20/455, 

wwd or wd6 s = 70/455, 

wwr or wrr t = 96/455, 

bbr ox brr u = i35/455» 

wbr V = 120/455. 

Prob. 7. When three dice are thrown together, what is the prob- 
ability that the throw will be greater than 9 ? 

Prob. 8. Write dcvn a literal formula for the probabilities of the 
several possible triads considered in the above question of the balls, 
supposing the numbers of white, black, and red balls to.be /, m, «, 
respectively. 
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Art. 5. Probability of Concurrent Events. 

If the probabilities of two independent events are/j and 
/a, respectively, the probability of their concurrence in any 
single instance is pj>^. Thus, suppose there are two urns 
f/j and f/j, the first of which contains ^j white and b^ black 
balls, and the second a^ white and b^ black balls. Then the 
probabiHty of drawing a white ball from f/j is/j = ajia^ -j" ^i)> 
while that of drawing a white ball from U^ is/, = aJ(a^-\- b^. 
The total number of different pairs of balls which can be formed 
from the entire number of balls is {a^ -\- b^{a^ + ^a)* O^ these 
pairs a^a^ are favorable to the concurrence of white in simul- 
taneous or successive drawings from the two urns. Hence the 
probability of a concurrence of 

white with white = ^a/(^, + ^iX^a + ^a)* 

white with black = {a^b^-^- aj?^/{a^ + *i)(^a + *a)> 

black with black = bJ)J{a^ + ^iX^a "h ^a)> 

and the sum of these is unity, as required by equations (2) of 
Article 4. 

In general, if A» A* A • • • denote the probabilities of several 
independent events, and P denote the probabiHty of theii 
concurrence, 

P^P.P.P.'*^ (i) 

To illustrate this formula, suppose there is required the 
probability of getting three aces with three dice thrown simul- 
taneously. In this case A =A =A = ^/^ ^^^ 
P={i/6y= 1/216. 

Similarly, if two dice arc thrown simultaneously the proba- 
bility that the sum of the numbers shown will be 1 1 is 2/36; 
and the probability that this sum ii will appear in two succes- 
sive throws of the same pair of dice is 4/36.36. 

The probability that the alternatives of a series of events 
will concur is evidently given by 

Q = q^q.q. . . . = (i - A)(i - A)(i -A) • • • (2) 

Thus, in the case of the three dice mentioned above, the 
probability that each will show something other than an ace is 
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q =i ^^ = /^^ =; 5^ 6, and the probaoiiicy "iiaX, rhe^/ -*iil concur in 

tiii> :-. (2 =^ i^5> -J^ 

.\rany cases of interest iccar for riie appiicarian of l .uid 
' 2K 0?ic of the most imoorrant jf riiesc :- riimisiied bv -uc* 
cessive tnais of the d&me event. C.jnsitien for tixarapie, '.viiiii 
niav hanoen in ft thais of an event for .viiicii :he ;inoaniiitv 
is p and a'^ainsC which the probaaiiir.' is /. T'le pnjbiibiiit\ 
that the e\'ent will occur everv rime ii y*. Tlie -irobaoiiitv lIiiiT 
the event .viil occur •/« — i) rimes in iu.cc;L5=ion and rhea f»al i*? 
/• ' v/. But if the order of occurre:ice is disre^^arded rhis liisr 
Ci>maination may arrive i\\ n iicFerent -.vays: -xk rhar rhe prob- 
abih'tv that the ^v*txiz will occur a — i rimes and fail once is 
up* - '//. Simiiariy, the pr:babiiiry tiiat the event will happei^ 
I y« — 2) rimes and faii twice is hn\ // — \) p ~ j \ etc That: is^ 
tiic orobabiiities of the several ccssibie occurrences are jven b'' 
tile corresponding terms in rhe development s:ii y^ — "if- 

By the same reasoning used to get equations kz\ oi Art 
I it may be shown thar the maximum term in the expansion 
of (p -^ qf^ is that in which the exponent m^ ssy, of ^ is 
t!ie whole number lying between w/-^ i)i^— t and l«-f- r)^* 
III other words, the mo^t probable result in « trials is the 
occurrence of the event (« — m) times and its failure m 
times. When /r is large this means that the most probable of 
all possible results is that in which the event occurs m ^ ftq 
z=. n(\ ^ q) ^ np times and fails «ytimesw Thus, if the event 
he that of throwing an ace with a single di^ the most probable 
of the p*->ssible results \n 6oo throws is that of icx> aces and 
y/j failures. 

Sinre q"* is the probability that the event will fail every time 
in n tri;il^ the probability that it will occur at least once in « 
trials is I — f. Calling this probability r* 

r- I -^= I -(!-/)-. ' (3> 

If r in this equation be replaced by I '2, the corresponding 

value of n is the number of trials essential to render the 

♦ Sk« Po!!W*'>n'4 Probability dcs Jogcmcnts, pp. ^o* 41. 
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chances even that the event whose probability is/ will occur 
at least once. Thus, in this case, the value of n is given by 

log(i-/)* 

This shows, for example, if the event be the throwing of double 
sixes with two dice, for which / = 1/36, that the chances are 
even {r = 1/2) that in 25 throws (« = 24.614 by the formula) 
double sixes will appear at least once. 

Equation (3) shows that however small/ may be, so long as 
it is finite, 71 may be taken so large as to make r approach in- 
definitely near to unity ; that is, n may be so large as to render 
it practically certain that the event will occur at least once. 

When n is large 

-^ ' 1.2 1.2.3' 
= e'"*^ approximately. 

Thus an approximate value of r is 

r = I - ^ - "^ log ^ = 0.4342495- (4) 

This formula gives, for example, for the probability of drawing 
the ace of spades from a pack of fifty-two cards at least once in 
104 trials r = I — ^"* = 0.865, while the exact formula (3) 
gives 0.867. 

Similarly, the probability of the occurrence of the event at 
least / times in n trials will be given by the sum of the terms 
of (/ -{- gY from p*" up to that in fq""^ inclusive. This proba- 
bility must be carefully distinguished from the probability that 
the event will occur / times only in the n trials, the latter being 
expressed by the single term in /'^""^ 

Prob. 9. Compare the probability of holding exactly four aces in 
five hands of whist with the probability of adding at least four aces 
in the same number of hands. 

Prob. 10. What is the probability of an event if the chances are 
even that it occurs at least once in a million trials? See equation (4). 
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t 

■ XK 

Art. 6. Bernoulli's Theorem. A 

Denote the exponents of/ and q in the maximum term of 
(/ + ^y ^y /* ^^^^ ^'^ respectively, and denote this term by 71 
Then 

As shown in Art. 5, /^ in this formula is the greatest whole 
number in (« + i)/, and m the greatest whole number in 
(« -[" i)^ ; so that when n. is large, fx and m are sensibly equal 
to np and ;/^ respectively. 

The direct calculation of 7 by (i) is impracticable when n 
is large. To overcome this difficulty the following expression 
is used : * 

nl = «V-'*i/2^/(i + — H--^5^ri+ . . .). (2) 

\ ' I2« ' 288« ' ' ^ ^ 

log e = 0.4342495, log 2;r = 0.7981799. 



This expression approaches //V"** S/inti as a limit with the 
y increase of n, and in this approximate form is known as Stir- 
ling's theorem. Although a rude approximation to ;/ ! for 
small values of n this theorem suffices in nearly all cases 
wherein such probabilities as 7 are desired. Making use of 

the theorem in (i) it becomes 

I 

^="-7 =• (3) 

ViTtnpq 

That this formula affords a fair approximation even when 
n is small is seen from the case of a die thrown 1/ times. The 
probability that any particular face will appear in one throw is 
/ = 1/6, whence q = 5/6; and the most probable result in 12 
throws is that .in which the particular face appears twice and 
fails to appear ten times. The probability of this result com- 
puted from (3) is 0.309, while the exact formula (i) gives 0.296. 

The probability that the event will occur a number of times 

* This expression is due to Laplace, Th^orie Analytique des Probabilit6s. 
See also De Morgan's Calculus, pp. 600-604. 
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« 

comprised between (/i — a) and (/i + ol) in « trials is evidently 
expressed by the sum of the terms in (/ + ?)" ^OJ* which the 
exponent of /has the specified range of values. Calling this 
probabiHty R^ putting 

fx =^ np -\- Uy and m = nq — «, 

and using Stirling's theorem (which implies that « is a large 
number),* 

R ^ :2—=====r\i -^ ■—] 1 1 

'\/2ntipq^ npi \ nql 

very nearly; and the summation is with respect to u from 
u ^=^ — a \.o u ^=^ -\- a. But expansion shows that the natural 
logarithm of the product of the two binomial factors in this 
equation is approximately — u^/2npq. Hence 

R := :^ ' ^-^^/^npq . 
^2nnpq 

and, since n is supposed large, this may be replaced by a definite 
integral, putting 



Thus 



dz = l/V2npq, and z^ = u^/2npq. 



-f o/ ^%npq a/ Vznpq 



^ 

9 



This equation expresses the theorem of James Bernoulli, 
given in his Ars Conjectandi, published in 171 3. 

The value of the right-hand member of (4) varies, as it 
should, between o and i, and approaches the latter limit rap- 
idly as s increases. Thus, writing for brevity 



z 



z 







* See Bertrand, Calcul des Probabilit^s, Paris, 1.889, for an extended discus- 
sion of the questions considered in this Article. 
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the following table shows the march of the integral : 



z 


/ 


z 


/ 


z 


/ 


O.OO 


O.CXX) 


0.75 


0.7II 


1.50 


o.q66 


.25 


.276 


1. 00 


.843 


1-75 


.987 


.50 


.520 


1.25 


.923 


2.00 


•995 






To illustrate the use of (4), suppose there is required the 
probability that in 6000 throws of a die the ace will ap{3ear a 
number' of times which shall be greater than 1/6 X 6000 — 10 
and less than 1/6 X 6000 -|- 10, or a number of times lying 
between 990 and loio. In this case a = 10, « = 6000,/ = 1/6, 

q =z 5/6. Thus, a/^2npq = 10/1^2 . 6000 . 1/6 . 5/6 = 0.245. - ^ r ,- ^ 

Hence, by (4) and the table, R = O.27. 

Prob. II. If the ratio of males to females at birth is 105 to loo," 
what is the probability that in the next 10,000 births the number of 
males will fall within two per cent of the most probable number ? 

Prob. 12. If the chance is even for head and tail in tossing a 
coin, what is the probability that in a million throws the difference 
between heads and tails will exceed 1500? 



?. 



Art. 7. Inverse Probabilities.* ^ 

If an observed event can be attributed to any one of several 
causes, what is the probability that any particular one of these 
causes produced the event ? To put the question in a concrete 
form, suppose a white ball has been drawn from one of two 
urns, U^ containing 3 white and 5 black balls, and C/, contain- 
ing 2 white and 4 black balls ; and that the probability in favor 
of each urn is required. If V^ is as likely to have been 
chosen as C/3,,the probability that C/, was chosen is 1/2. After 
such choice the probability of drawing a white ball from C/, is 
3/8. Before drawing, therefore, the probability of getting a 
white ball from t/j was 1/2 X 3/8 = 3/16, by Art. 5. Similarly, 
before drawing the probability of getting a white ball from U^ 
was 1/2 X 2/6 = 1/6. These probabilities will remain un- 
changed if the number of balls in either urn be increased or 

* See Poisson, Probability des Jugements, pp. 81-83. 
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diminished so long as the ratio of white to black balls is kept 
constant. Make these numbers the same for the two urns. 
Thus let the first contain 9 white and 15 black, and the second 
8 white and 16 black ; whence the above probabilities may be 
written 1/2 x 9/24 and 1/2 X 8/24. It is now seen that there 
are (9 -|- 8) cases favorable to the production of a white ball, 
each of which has the same antecedent probability, namely, 1/2. 
Since the fact that a white ball was drawn excludes considera- 
tion of the black balls, the probability that the white ball came 
from U^ is 9/17 and that it came from U^ is 8/17; and the sum 
of these is unity, as it should be. 

To generalize this result, let there be m calises, C,, C„ . . . C^. 
Denote their direct probabilities by ^,,^a,. ..^„,; their antecedent ;< 
probabilities by r,, r^, . . .r^\ and their resultant probabilities 
on the supposition of separate existence by p^^p^y . . ./,»• 
That is, 

Px = ?i^» A = ?a^» ' * 'Pm = ^m^m- (0 

Let Z? be the common denominator of the right-hand mem- 
bers in (i), and denote the corresponding numerators of the 
several fractions by s.,s^,.., s,^. Then 

p, = sJD, p^ = sJD, .../,«= sJD ; 

and it is seen that there are in all (^, + ^, + . . . s^ equally 
possible cases, and that of these s^ are favorable to C,, j, to 
C,, . . . Hence, if /\, /^i, . . . Pnt denote the probabilities of , 
the several causes on the supposition of their coexistence, 

Py = ^i/(^. + ^. + • • • -^m) ^PAP, + a + . . . /J. 

Thus in general 

P. = /,/^A P. = A/^A .../"«.= A/^A (2) 

To illustrate the meaning of these formulas by the above 
concrete case of the urns it suffices to observe that 

for 6^„ q^ = 3/8 and r, = 1/2, 

for C/,, ^, = 1/3 and r^ = 1/2 ; 
whence /, = 3/16, p, = 1/6. /^ +/, = 17/48 ; 

and /», =9/17, P^ = S/\7. 

As a second illustration, suppose it is known that a white 



26 PROBABILITY AND THEORY OF ERRORS. 

ball has been drawn from an urn which originally contained in 
balls, some of them being black, if all are not white. What is 
the probability that the urn contained exactly n white balls? 
The facts are consistent with in different and equally probable 
hypotheses (or causes), namely, that there were i white and 
(ill — i) black balls, 2 white and (/// — 2) black balls, etc. 
Hence in (i), ^^ = ^^ = . . . = i, and 

p^ = i/nty /a = 2/w, .../»„ = n/m, . . . /^ = in/m. 

Thus 2p = (i/2){m + i), 

2;/ 
and Pn ^=^pn/2p = ~~, ; — :. 

This shows, as it evidently should, that n = in is the most 
probable number of white balls in the urn. The probability 
for this number is P^ = 2/(m -\- i), which reduces, as it ought, 
to I for in= I, 

Formulas (i) and (2) may also be applied to the problem of 
estimating the probability of the occurrence of an event from 
the concurrent testimony of several witnesses, Jf,, Jf,, . . . 
Denote the probabilities that the witnesses tell the truth by 
x^ x^y , . , Then, supposing them to testify independently, 
the probability that they will concur in the truth concerning 
the event is x^x^ . . .; while the probability that they will con- 
cur in the only other alternative, falsehood, is (i —x^){i —x^ . . . 
The two alternatives are equally possible. Hence by equations 
(i) and (2) 

/. = ^.-^^ • • M A = (i - -^iXi - -^a) • • •» 

p -^i-^a * * • 

* ~ x,x^ . . . -f (l — x,){\ - .ar,) . . .' 

(l-x ,){l^ x:)... ^3) 

« X,X^ . . . + (l — X^)( I — x^) . . .* 

P^ being the probability for and P^ that against the event. 

To illustrate (3), if the chances are 3 to I that X^ tells the 
truth and 5 to i that X, tells the truth, x^ = 3/4, x^ = 5/6, and 
Pi= 15/16; or, the chances are 1$ to i that an event occurred 
if they agree in asserting that it did.* 

* For some interesting applications of equations (3) see note E of Appendix 
to the Ninth Bridgewater Treatise by Charles Babbage (London, 1838). 
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It is of theoretical interest to observe that if x^t x^, . . , \n 

(3) are each greater than 1/2, /\ approaches unity as the 

number of witnesses is indefinitely increased. 

Prob. 13. The groups of numbers of one figure each, two figures 
each, three figures each, etc., which it is possible to form from the 
nine digits i, 2, ... 9 are printed on cards and placed severally in 
nine similar urns. What is the probability that the number 777 will 
be drawn in a single trial by a person unaware of the contents of 
the urns ? 

Prob. 14. How many witnesses whose credibilities are each 3/4 ^, 
are essential to make P^ = 0.999 ^^ equation (3) ? 

Art. 8. Probabilities of Future Events. 

Equations (2) of Art. 7 may be written in the following 
manner: „ ,^ „ . 

A A A ^/ ^'^ 

If/,, A» • • • P»t are found by observation, P^, /!,,... P„t will ex- 
press the probabilities of the corresponding causes or their 
effects. When, as in the case of most physical facts, the num- 
ber of causes and events is indefinitely great, the value of any 
p or P in (i) becomes indefinitely small, and the value of 2p 
must be expressed by means of a definite integral. Let x de- 
note the probability of any particular cause, or of the event to - 
which it gives rise. Then, supposing this and all the other 
causes mutually exclusive, (i — x) will be the probability 
against the event. Now suppose it has been observed that in 
(;// -|- n) cases the event in question has occurred in times and 
failed n times. The probability of such a concurrence is, by 
Art. 5, cx\\ — xY, where c is a constant. Since x is unknown, 
it may be assumed to have any value within the limits o and I ; 
and all such values are k priori equally possible. Put 

y = cx^^'^i — xy. 

Then evidently the probability that x will fall within any as- 
signed possible limits a and h is expressed by the fraction 



O 1 

Jydx I j ydx\ 
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SO that the probability of any particular x is given by 



Cx\\ — x^dx 



This may be regarded as the antecedent probability of the 
cause or event in question. 

What then is the probability that in the next (r + s) trials 
the event will occur r times and fail s times, if no regard is had 
of the order of occurrence? If x were known, the answer 
would be by Arts. 2 and 5 

-Sr^^('-^>'- (3) 

But since x is restricted only by the condition (2), the required 
probability will be found by taking the product of (2) and (3) 
and integrating throughout the range of x. Thus, calling the 
required probability Q, 

^ r\s\ I • ^4> 

/ x\\ — xydx 



The definite integrals which appear here are known as Gamma 
functions. They are discussed in all of the higher treatises on 
the Integral Calculus. Applying the rules derived in such 
treatises there results * 

^"" r\s\m\7t\{fn^n-\-r'Ys'\-\)\ * ^5} 

If regard is had to the order of occurrence of the event ; 
that is, if the probability required is that of the event happen- 
ing r times in succession and then failing s times in succession, 

* It is a remarkable fact that formula (5) is true without restriction as to 
values of w, n, r, j. The formula may be established by elementary considera- 
tions, as was done by Prevost and Lhuilier, 1795. See Todhunter's History of 
he Theory of Probability, pp. 453-457- 
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the factor {r-^s)l/r\sl in (3), (4), (5) must be replaced by 
unity. 

To illustrate these formulas, suppose first that the event 
has happened m times and failed no times. What is the prol>- 
ability that it will occur at the next trial? In this case (4) 
gives 

Q = fx'^'^'dx I Jx'^dx = (w + \)/{fn + 2). 



When m is large this probability is nearly unity. Thus, the 
sun has risen without failure a great number of times m ; the 
probability that it will rise to-morrow is 

\ ^ mf\ ^ ml m ^ itr 

which is practically i. 

Secondly, suppose an urn contains white and black balls in 
an unknown ratio. If in ten trials 7 white and 3 black balls 
are drawn, what is the probability that in the next five trials 
2 white and 3 black balls will be drawn ? The application of 
(5) supposes the ratio of the white and black balls in the urn 
to remain constant. This will follow if the balls are replaced 
after each drawing, or if the number of balls in the urn is sup- 
posed infinite. The data give 

w = 7, « = 3» r = 2, ^ = 3» 

m-\-7i-\-r-\-s-\'\ = 16. 
Thus by (5) 

n 5!9!6!ii! , 

^ =2!3!7!3!i6! = ^^/^l 

Suppose there are two mutually exclusive events, the first 
of which has happened m times and the second n times in 
m -\-n trials. What is the probability that the chance of the 
occurrence of the first exceeds 1/2 ? The answer to this ques- 
tion is given directly by equation (2) by integrating the nume- 
rator between the specified limits of x. That is. 



< 



>- 
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y ;r'"(l — xfdx 

P^-. . (6) 

I x*^{i —xfdx 

• 
Thus, it m = I and k ^o, P = 3/4 ; or the odds are three to 
one that the event is more likely to happen than not. Simi- 
larly, if the event has occurred w times in succession, 

which approaches unity rapidly with increase of n. 

Art. 9. Theory of Errors. 

The theory of errors may be defined as that branch of math- 
ematics which is concerned, first, with the expression of the re- 
sultant effect of one or more sources of error to which com- 
puted and observed quantities are subject ; and, secondly, with 
the determination of the relation between the magnitude of 
an e^rror and the probability of its occurrence. In the case of 
computed quantities which depend on numerical data, such as 
tables of logarithms, trigonometric functions, etc., it is usually 
possible to ascertain the actual values of the resultant errors. 
In the case of observed quantities, on the other hand, it is not 
generally possible to evaluate the resultant actifal error, since 
the actual errors of observation are usually unknown. In either 
case, however, it is always possible to write down a symbolical 
expression which will show how different sources of error enter 
and affect the aggregate error; and the statement of such an 
expression is of fundamental importance in the theory of errors. 

1 To fix the ideas, suppose a quantity Q to he a function of 
several independent quantities ;ir, j, ^ . . . ; that is, 

and let it be required to determine the error in Q due to errors 
in X, j^f js . . . Denote such errors by JQ, Ax, Ay^ Az . , . 
Then, supposing the errors so small that their squares, prod- 
ucts, and higher powers may be neglected, Taylor's series gives 
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.^ dQ . , dQ. , dQ^ , 
^G=9^^^ + ^^J'+^^^ + ... (I) 

This equation may be said to express the resultant actual error 
of the function in terms of the component actual errors, since 
the actual value of JQ is known when the actual errors of 
x,y, z , . . are known. It should be carefully noted that the 
quantities x^y^ z . .. are supposed subject to errors which are 
independent of one another. The discovery of the independent 
sources of error is sometimes a matter of difficulty, and in general 
requires close attention on the part of the student if he would 
avoid blunders and misconceptions. Every investigator in work 
of precision should have a clear notion of the error-equation of 
the type (i) appertaining to his work; for it is thus only that 
he can distinguish between the important and unimportant 
sources of error. 

Prob. 15. Write out the error-equation in accordance with (i) 
for the function Q -= xyz + x^ log {y/z), 

Prob. 16. In a plane triangle a/b = sin A /sin B, Find the error 
in a due to errors in ^, //, and B. 

Prob. 17. Suppose in place of the data of problem i6 that the 
angles used in computation are given by the following equations : 
^=^.+i(i8o°-^,- ^.- C,), B=B,+ i(i8o° - A, -^. - CX 
where A^, B^, C, are observed values. What then is Aa? 

Prob. 18. If 7a denote the weight of a body and r the radius of 
the earth, show that for small changes in altitude, Aw/w^=- — 2Ar/r\ 
whence, if a precision of one part in 500 000 000 is attainable in com- 
paring two nearly equal masses, the effect of a difference in altitude 
of one centimeter in the scale-pans of a balance will be noticeable.* 

Art. 10. Laws of Error. 

A law of error is a function which expresses the relative 
frequency of occurrence of errors in terms of their magnitudes. 
Thus, using the customary notation, let e denote the magni- 

* This problem arose with the International Bureau of Weif? hts and Measures, 
whose work of intercomparison of the Prototype Kilogrammes attained a pre- 
cision indicated by a probable error of 1/500 000 000th part of a kilogramme. 
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tude of any error in a system of possible errors. Then the law 
of such system may be expressed by an equation of the form 

y = 0(6). (I) 

Representing e as abscissa and y as ordinate, this equation 
gives a curve called the curve of frequency, the nature of which, 
as is evident, depends on the form of the function 0. This 
equation gives the relative frequency of occurrence of errors in 
the system ; so that if e is continuous the probability of the 
occurrence of any particular error is expressed hy yds = 0(e)^/e; 
which is infinitesimal, as it plainly should be, since in any con- 
tinuous system the number of different values of e is infinite. 

Consider the simplest form of 0(e), namely, that in which 
0(e) = ^, a constant. This form of 0(6) obtains in the case of 
the errors of tabular logarithms, natural trigonometric func- 
tions, etc. In this case all errors between minus a half-unit 
and plus a half-unit of the last tabular place are equally likely 
to occur. Suppose, to cover the class of cases to which that 
just cited belongs, all errors between the limits — ^ and -\-a 
are equally likely to occur. The probability of any individual 
error will then be 0(e)^/e =- cde, and the sum of all such prob- 
abilities, by equation (2), Art. 4, must be unity. That is, 

I (p{€)d€ = c ide = I. (2) 



— a 



—a 



This gives c = 1/2^, or by (i) y = \/2a. The curve of fre- 
quency in this case is shown in the figure, 
AB being the axis of e and OQ that of ^. 
It is evident from this diagram that if the 
errors of the system be considered with 
respect to magnitude only, half of them 
should be greater and half less than a/2. 
This is easily found to be so in the case of 
tabular logarithms, etc. 

As a second illustration of (i), suppose j/ and e connected 

by the relation j^ = c \/a* — e', where a is the radius of a circle, 



c, j> 

A o B 
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c a constant, and e may have any value between — a and + ^ 
Then the condition 

cjde )/'^^^' = I 



- a 



gives c = 2/{a^7i), In this, as in the preceding case, 0(-|- e) = 
<p{— e), the meaning of which is that positive and negative- 
errors of the same magnitude are equally likely to occur. It 
will be noticed, however, that in the latter case small errors 
have a much higher probability than those near the limit a, 
while in the former case all errors have the same probability. 

In general, when e is continuous 0(e) must satisfy the condi- 
tion / (t)(e)de = I, the limits being such as to cover the entire 

range of values of e. The cases most commonly met with are 
those in which 0(6) is an even function, or those in which 
0(-|- e) = 0(— e). In such cases, if ± ^ denote the limiting 
value of €, 

-fa a 

J'(p{e)de = 2j<p{e)cle = I. (3) 



— a 



Art. 11. Typical Errors of a System. 

Certain typical errors of a system have received special 
designations and are of constant use in the literature of the 
theory of errors. These special errors are the probable error, 
the mean error, and the average error. The first is that error 
of the system of errors which is as likely to be exceeded as 
not ; the second is the square root of the mean of the squares 
of all the errors ; and the third is the mean of all the errors 
regardless of their signs. Confining attention to systems in 
which positive and negative errors of the same magnitude 
are equally probable, these typical errors are defined mathe- 
matically as follows. Let 

€p = the probable error, 

€„,= the mean error, 

Ea = the average error. 
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Then, observing (2j, of Art. lo, 

f<He)de =f<p{e)de = J <p(e)de =f4>(e)de = f 

€^" = / <ti{e)^de^ e^ = 2 / 0(e) e^e. 



— a 



(0 



The student should seek to avoid the very common misap- 
prehension of the meaning of the probable error. It is not 
** the most probable error," nor ** the most probable value of 
the actual error"; but it is that error which, disregarding signs, 
would occupy the middle place if all the errors of the system 
were arranged in order of magnitude. A few illustrations will 
suffice to fix the ideas as to the typical errors. Thus, take the 
simple case wherein 0(e) = ^ = 1/2^, which applies to tabular 
logarithms, etc. Equations (i) give at once 

e^ = ± -a, €^= ± -y i, ea= ±-a. 
2 3 2 

For the case of tabular values, a = 0.5 in units of the last 

tabular place. Hence for such values 

€p = ±0.25, e^ = ± 0.29, €^ - ± 0.25. 

Prob. 19. Find the typical errors for the cases in which the law 

of error is 0(e) = cVa^ — e*, 0(e) =c(±a^^e), (pie) =c cos* (Tie/ 2a); 
c being a constant to be determined in each case and e having any 
value between — a and + a. 

Art. 12. Laws of Resultant Error. 

When several independent sources of error conspire to pro- 
duce a resultant error, as specified by equation (i) of Art. 9, 
there is presented the problem of determining the law of the 
resultant' error by means of the laws of the component errors. 
The algebraic statement of this problem is obtained as follows 
for the case of continuous errors : 

In the equation (i). Art. 9, write for brevity 

0^ dy 
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and let the laws of error of e, e^, e,, . . . be denoted by 0(e), 
j^j(ej, 0a(e,) . . . Then the value of e is given by 

e=e, + e,+. .. (l) 

The probabilities of the occurrence of any particular values 
of Cj, e„ . . . are given by 0j(e,)^e,, 0a(e,)^e„ . . . ; and the 
probability of their concurrence is the probability of the cor- 
responding value of €, But since this value may arise in an 
infinite number of ways through the variations of e,, e„ . . . 
over their ranges, the probabiHty of e, or 0(e)^e, will be 
expressed by the integral of (pi(€i)de^(p^{€ ^de^ . . . subject to 
the restriction (i). This latter gives e, = e— e^ — e, . . ., and 
de^ = de for the multiple integration with respect to e,, e„ . . . 
Hence there results 

(p{e)de = de C (t>S^e — e, — e, — . . .)0a(^a)^*^a • • • > 



or 



0(^) - f<t>X^-^x - ^a - • • ')<t>J,eye^f(t),{e,)de^ ... (2) 

It is readily seen that this formula will increase rapidly in 
complexity with the number of independent sources of error.* 
For some of the most important practical applications, how- 
ever, it suffices to limit equation (2) to the case of two inde- 
pendent sources of error, each of constant probability within 
assigned limits. Thus, to consider this case, let e, vary over 
the range — a to -\-ay and e, vary over the range — ^ to + ^« 
Then by equation (2), Art. 10, 

0,(€,) = 1/(2^), 0,(6,) = 1/(2*). 

Hence equation (2) becomes 

In evaluating this integral e, must not surpass ± ^ and 
e, = e — €3 must not surpass ± ^^ Assuming a> 6y the limits 
of the integral for any value of e = e, + e, lying between 

— (^ + ^) and —{a — 6) are — d and + (e -|- a). This fact is 

I 

* The reader desirous of pursuing this phase of the subject should consult 
Bessel's Uritersuchungen ueber die Wahrscheinlichkeit der Beobachtungsfehler; 
Abhandlungen von Bessel (Leipzig, 1876), VoL IL 
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made plain by a numerical example. For instance, suppose 
a = $ and ^ = 3. Then — (^^ + ^) = — 8 and -"{a— 6)= — 3. 
Take e = — 6, a number intermediate to — 8 and — 3. Then 
the following are the possible integer values of c, and c, which 
will produce 6 = — 6 : 

€ €, €, limits of €, 

-6^-5-1, -!= + (€ + ^jX 
= - 4 - 2, 

= -3-3, -3=-*- 

Similarly, the limits of c, for values of e lying between 
—■ {a — b) and + (^ — ^) are — b and + ^ ; and the limits of 
€, for values of e between -f (a — ^) and -|- {a-\- b) are -}-(€--/?) 
and + b. Hence 

•+« 



^^)=^^ A-= ji for - (— ^ <-<+(— *). 



-* 
+* 



^^>=^*/''^'=''^' ^^-^ +(-*)<^<+(-+^)- 



(3) 



•-« 



Thus it appears that in this case the graph of the resultant 
law of error is represented by the upper base and the two sides 

of a trapezoid, the lower base beinor the 

IIOII f » i> 

axis of 6 and the line joining the middle 

points of the bases being the axis of 0(e). 
Ill loi III 

(See the first figure in Art. 13.) The prop- 
in IIOI IIII r-r- 

erties of (3), including the determination 

IIIIIIOIIIIII , , ,. . 1 11 . . 1 1 

of the limits, are also illustrated by the 
IIIIIIIOIIIIIII ,. , , ., f 1 , 

adiacent trapezoid of numerals arranged 
IIIIIIIIOIIIIIIII •' /. , , ^ , 

to represent the case wherein a = 0.5 and 
^ = 0.3. The vertical scale, or that for 0(e), does not, how- 
ever, conform exactly to that for e. 
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Prob. 20. Prove that the values of 0(6) as given by equation (3) 
satisfy the condition specified in equation {3), Art. 10. 

Prob. 21. Examine equations {3) for the cases wherein ^jr = ^and 
^ = o; and interpret for the latter case the first and last of (3). 

Prob. 22. Find from (3), and (i) of Art. 11, the probable error of 
the sum of two tabular logarithms. 

Art. 13. Errors of Interpolated Values. 

Case I. — One of the most instructive cases to which formulas 
(3) of Art. 12 are applicable is that of interpolated logarithms, 
trigonometric functions, etc., dependent on first differences. 
Thus, suppose that v^ and v^ are two tabular logarithms, and 
that it is required to get a value v lying / tenths of the interval 
from r, towards v^. Evidently 

^ = ^i + (^a — ^1) ^ = (i — ^M + i'^1 ; 

and hence if ^, e^, e^ denote the actual errors of v, v^, v^ , re- 
spectively, 

e = {i-t)e, + te,. (i) 

It is to be carefully noted here that e as given by (i) re- 
quires the retention in v of at least one decimal place be- 
yond the last tabular place. For example, let v = log (24373) 
from a 5-place table. Then t\ = 4.38686, 7/, = 4.38703, 
v^ — z/, = +0.00017, / = 0.3, and 7f =4.38691.1. Likewise, as 
found from a 7-place table, e^ •= —0.45, e.^ = -[-0.37 in units of 
the fifth place; and hence by (i) e= —0.20. That is, the 
actual error of z/ = 4.38691. 1 is = 0.20, and this is verified by 
reference to a 7-place table. 

The reader is also cautioned against mistaking the species 
of interpolated values here considered for the species common- 
ly used by computers, namely, that in which the interpolated 
value is rounded to the nearest unit of the last tabular place. 
The latter species is discussed under Case II below. 

Confining attention now to the class of errors specified by 

equation (i), there result in the notation of the preceding 

article 

6, = (l— /)^p 6a = /^a, and € = ^=6, + e, ; 

and since e^ and e^ each vary continuously between the limits 
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for values of e betw. —0.5 and —(0.5—/), 



± 0.5 of a unit of the last tabular place, a and b in equations 
(3) of that article have the values 

a = 0.5(1 — /), b = 0.5/. 

Hence the law of error of the interpolated values is ex- 
pressed as follows : 

0.5 + e 

= for values of e betw. —(0.5—/) and +(0.5—/), \ (2) 

= -^ for values of e betw. +(0.5 — /) and +0.^. 

The graph of 0(e) for / = 1/3 is shown by the trapezoid 
ABf BC^ CD in the figure on page 40. Evidently the equa- 
tions (2) are in general represented by a trapezoid, which degen- 
erates to an isosceles triangle when / = 1/2. 

The probable, mean, and average errors of an interpolated 
value of the kind in question are readily found from (2), and 
from equations (i) of Art. 11, to be 

ep = (i/4)(i - /) for o < / < 1/3. ' 

= 1/2 - (1/2)1^2/(1 -/) 
= 1/4/ 



for 

for 1/3 < / < 2/3, 
for 2/3 </ < 1. 



e- = 



'W 



Ca = 



_ ( i -d -2/)' )'/' 

- I 96(1 - /y ) * 

I - (I - 2/)' 



y (3) 



24(1 — /)/ 

I _ (2/ -^ I)' 

24(1 — /)/ 



for o < / < 1/2, 
for 1/2 </ < I. 



It is thus seen that the probable error of the interpolated 
value here considered decreases from c.25 to 0.15 of a unit of 
the last tabular place as / increases from o to 0.5. Hence such 
values are more precise than tabular values; and the computer 
who desires to secure the highest attainable precision with a 
given table of logarithms should retain ^ne additional figure 
beyond the last tabular place in interpolated values. 
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Case 11.^ — Recurring to the equation 2/ = v^-\- i{v^ — v^) for an 
interpolated value v in terms of two consecutive tabular values 
v^ and z/j, it will be observed that if the quantity /(z/, — v^) is 
rounded to the nearest unit of the last tabular place, a new error 
is introduced.. For example, if z/j = log 1633 = 3.21299, and 
v^ = log 1634 = 3.21325 from a 5-place table, v, — v^= -{- 26 
units of the last tabular place ; and if / = 1/3, /(z;, — v^) = 8| ; 
so that by the method of interpolation in question there results 
V = 3.21299 +9= 3.21308. Now the actual errors of t\ and 
v^ are, as found froma7-place table, — 0.38 and -|-0 2i in units 
of the fifth place. Hence the actual error of v is by equation 
(0> f X — 0.38 + i X + 0.21 — J = — 0.52, as is shown di- 
rectly by a 7-place table. \ 

It appears, then, that in this case the error-equati6n cor- 
responding to (i) is 

^= (i-/>, + /^, +^., \ (4) 

wherein e^ and ^, are the same as in (i) and ^, is the actualerror 
that comes from rounding /(v^ — v^) to tlie nearest unit oif the 
last tabular place. 

The error e^, however, differs radically in kind from e^ and 
e^. The two latter are continuous, that is, they rnay each have 
any value, between the limits — 0.5 and +0-5 J while ^, is dis- 
continuous, being limited to a finite number of values depend- 
ent on the interpolating factor /. Thus, for / = 1/2 the only 
possible values of ^, are o -(- 1/2, and — 1/2 ; likewise for / =. 
1/3, the only possible values of e^ are o, + i/3» ^"d — i/S* ^^ 
is also clear that the maximum value of e, which is constant and 
equal to 1/2 for (i), is variable for (4) in a manner dependent 
on /. For example, in (4), 

The maximum o{ e = 1/2 + 1/2 = i, for / = 1/2, 

" ^ = 1/2 + 1/3 = 5/6, " / =r 1/3, 
" ^ = 1/2 + 1/2 =1, " / = 1/4, 
" ^ = 1/2 + 2/5 =9/10 " / =: 1/5. 

The determination of the law of error for this case presents 
some novelty, since it is essential to combii^e the continuous 
errors (i ^ /y^ and fe^ with the discontinuous error e^. The 
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simplest mode of attacking the problem seems to be the fol- 
lowing quasi-geometrical one. In the notation of Arts. 12 and 
13, put in (4) e = e,{l — /V, := e,, te, = e,, and <■, = e,. Then 
e=K + e,) + e.- (5) 

The law of error for («, -|- €,) is given by equation (2) for any 
value of /. Hence for a given value of / there will be as many 
expressions of 0(e) as there are different values of e,. The 
graphs of 0(6) will at! be of the same form but will be differently 
placed with reference to the axis of 0(e). Thus, if / = 1/3 the 
values of «, are — 1/3, o, and 
+ i/3> ^"tl these are equally 
likely to occur. For e, = O the 
graph is given directly by (2), 
and is the trapezoid ABCD 
symmetrical with respect to OQ. 
For e, = — 1/3 the graph is 
abQd, of the same form as 
ABCD but shifted to the left 
by the amount of c, = — 1/3, 
Similarly, the graph for the case 
of «, = -f 1/3 is t^Qb'd', and is produced by shifting ABCD to 
the right by an amount equal to + 1/3. 

Now, since the three systems of errors for this case are 
equally likely to occur, they may be combined into one system 
by simple addition of the corresponding element areas of the 
several graphs. Inspection of the diagram shows* that the 
resultant law of error is expressed by 

<^«) = (f/4)(5 + 6e) for -5/6<6<- ,/6,^ 

= 1 for — 1/6 < e < + 1/6, \ (6) 

= (1/4XS - 6e) for + 1/6 < e < + 5/6. J 
This is represented by a trapezoid whose lower base is 10/6, 
upper base 2/6, and altitude 1. 




* Slim (he three areas and 
required by cquadon (3), Art. 



divide by 3 lo make resuliant ei 
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As a second illustration, consider equation (5) for the case 
/ = 1/2. In this case e, must be either o or 1/2, the sign of 
which latter is arbitrary. For e, = o, equations (2) give 

0(e) = 2 + 4€ for - 1/2 < e < o, I 

= 2 — 4e for o < e <+ 1/2. ) ^7) 

This function is represented by the isosceles triangle A QE 

whose altitude OQ is twice the base AE, 

Similarly 0(e) for e, = -j- 1/2 would 

be represented by the triangle -^4 (2-£ dis- 
placed to the right a distance 1/2 ; and 
if the two systems for e, = o and e, = 
-)- 1/2 be combined into one system, 
their resultant law of error is evidently 
0(e)= i-|-2e for — i/2<e< o, \ 

= I foro< e<+i/2, (. (8) 

= 2 — 2e for + 1/2 < e< I ; ) 
the graph of which is ABCD. On the 

other hand, if the errors in this combined system be considered 
with respect to magnitude only, the law of error is 

0(e) = 2(1 -e) for o < e< i, (9) 

the graph of which is OQD. 

The student should observe that (6), (7), (8), and (9) satisfy 

the condition / 0(e)<ie = i if the integration embraces the 
whole range of e. 

The determination of the general form of 0(e) in terms of 
the interpolating factor / for the present case presents some 
difficulties, and there does not appear to be any published solu- 
tion of this problem.* The results arising from one phase of 
the problem have been given, however, by the author in the 
Annals of Mathematics,! and may be here stated without 
proof. The phase in question is that wherein t is of the form 
i/«, n being any positive integer less than twice the greatest 

* The author explained a general method of solution in a paper read at the 
summer meeting of the American Mathematical Society, August, 1895. 
t Vol. II, pp. 54-59- 
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tabular difference of the table to which the formulas are ap- 
plied. For this restricted form of t the possible maximum 
value of e as given by equation (5) is, in units of the last 
tabular place, (2« — \)/n for n odd and i for n even. 

The possible values of e, of equation (5) are 

I 2 n— \ 

o, ± -, ± -, . . . ± for n odd, 

I 2 n — 2 I ^ 

o, ± — , ± -, . . . ± , ± - for ;^ even. 

/ » n 2/12 

An important fact with regard to the error 1/2 for n even 
is that its sign is arbitrary, or is not fixed by the computation 
as is the case with all the other errors. However, the com- 
puter's rule, which makes the rounded last figure of an inter- 
polated value even when half a unit is to be disposed of, will, 
in the long-run, make this error as often plus as minus. 

The laws of error which result are then as follows : 

For n odd. 

0(e) = I for € between — i/2n and + i/2^» 

0(e) = ( ± ej for ebetw. T i/2« and T(2«— 1)/2«. 

For n even. 

n /2n — 2 \ 

0(e) = —. A + el for e between o and T i/«, 

^ ^ 2(« — i)\ « / ' / ' . 

n (2n — I \ 

= 1 ± el for € betw. =F i/« and =F (« — i)/«» 

ti 

= — r(i ± e) for e between T («— lA) and T I. 

2{n — ly ^ 

By means of these formulas and (i) of Art. 1 1 the probable, 
mean, and average errors for any value of ;/ can be readily 
found. The following table contains the results of such a com- 
putation for values of n ranging from i to 10. The maximum 
actual error for each value of n is also added. The verifica- 
tion of the tabular quantities will afford a useful exercise to the 
student. 
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Typical Errors of Interpolated Logarithms, etc. 



Interpolating 
Factor. 


Probable Error. 


Mean Error. 


Average Error. 


Maximum 
Actual Error. 


I 


0.250 


0.289 


0.250 


:/2 


1/2 


.292 


.408 


.333 


I 


1/3 


.256 


.347 


.287 


S/6 


1/4 


.276 


.382 


.3>3 


I 


«/5 


.268 


•370 


•303 


9/10 


1/6 


277 


.385 


.315 


I 


1/7 


.274 


.380 


•3" 


13/14 


l/i 


.279 


.389 


.318 


I 


1/9 


.278 


.386 


.316 


17/18 


l/lO 


.281 


.392 


.320 


I 



When the interpolating factor / has the more general form 
m/n, wherein m and n are integers with no common factor, the 
possible values of e, are the same as for t = i/«. But equa- 
tions (3) of Art. 12 are not the same for / = m/^ as for/ = i/«, 
and hence for the more general form of t, 0(e) assumes a new 
type which is somewhat more complex than that discussed 
above. The limits of this work render it impossible to extend 
the investigation to these more complex forms of 0(e). It may 
suffice, therefore, to give a single instance of such a function, 
namely, that for which / = 2/5. For this case 

0(e) =1 for e between o and T i/io, 

= (5/6)(i3/iO ± e) for e between =F i/io and ip 3/10, 
= (5/3X4/5 ± e) for e between :f 3/10 and T 7/10, 
= (S/6)(9/io ± e) for e between T 7/iO and ^ 9/10. 
The graph of the right-hand half of ^ b 
this function is shown in the accompany- 
ing diagram, the whole graph being 
symmetrical with respect to OA, or the 
axis of 0(e). 

i. Attention maybe called to the strik- 
ing resemblance of this graph to that of 
the law of error of least squares. 

Prob. 23. Show from equations (3) that e^ varies from 1//12 
= 0.29 — , for / = o, to i/i^24 = 0.20 +, for / = 0.5 ; and that e« 
varies from 0.25 to 1/6 for the same limits. 
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Prob. 24. Show that the probable, mean, and average errors 
for the case of /= 2/5 cited above (p. 43) are ± 0.261, ± 0.251, 
and ± 0.290, respectively. 

Art. 14. Statistical Test of Theory. 

A statistical test of the theory developed in Art. 13 may 
be readily drawn from any considerable number of actual er- 
rors of interpolated values dependent on the same interpolating 
factor. The application of such a test, if carried out fully by 
the student, will go far also towards fixing clear notions as to 
the meaning of the critical errors. 

Consider first the case in which an interpolated value falls 
midway between two consecutive values, and suppose this 
interpolated value retains two additional figures beyond the 
last tabular place. Then by equations (2), Art. 13, the law of 
error of this interpolated value is 

0(e) = 2 -|- 46 for € between — 0.5 and o 
= 2 — 46 for 6 between o and + 0.5. 

Hence by equation (i) of Art. 1 1, or equation (3) of Art. 12, the 
probable error in this system of errors is i — (J) V2^= 0.15. 
It follows, therefore, that in any large number of actual errors 
of this system, half should be less and half greater than 0.15. 
Similarly, of the whole number of such errors the percentage 
falling between the values 0.0 and 0.2 should be 

J (l){€)de = 2y (2 — ^e)de = 0.64 ; 



-o.a 



that IS, sixty- four per cent of the errors in question should be 
less numerically than 0.2. 

To afford a more detailed comparison in this case, the act- 
ual errors of five hundred interpolated values from a 5-place 
table have been computed by means of a 7-place table. The 
arguments used were the following numbers : 20005, 20035, 
20065, 20105, 20135, etc., in the same order to 36635. The 
actual and theoretical percentages of the whole number of 
errors falling between the limits 0.0 and o.i, o.l and 0.2, etc., 
are shown in the tabular form following : 
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Limits of Errors. ^ ^^^"f ^ Theoretical 

Percentage. Percentage. 

o.o and o.i 33.2 36 

o. I and 0.2 30.2 28 

0.2 and 0.3 19.0 20 

0.3 and 0.4 13.2 12 

0.4 and 0.5 4.4 4 

0.0 and 0.1 5 51.4 50 

The agreement shown here between the actual and theoretical 
percentages is quite close, the maximum discrepancy being 2.8 
and the average 1.5 per cent. 

Secondly, consider the case of interpolated mid-values of the 
species treated under Case II of Art. 13, The law of error for 
this case is given by the single equation (9) of Art. 13, namely, 
0(6) = 2(1 — e), no regard being paid to the signs of the errors. 
The probable error is then found from 



• 



whence e^ = i — i i^ = 0.29. Similarly, the percentage of 
the whole number of errors which may be expected to lie, for 
example, between 0.0 and 0.2 in this system is 



9 



2 / (l — €)d€ = 0.36. 



• 



Using the same five hundred interpolated values cited 
above, but rounding them to the nearest unit of the last tabu- 
lar place and computing their actual errors by means of a 7-place 
table, the following comparison is afforded : 

T . .. t IT „ Actual Theoretical 

Limits of Errors. Percentage. Percentage. 

0.0 and 0.2 35.8 36 

0.2 and 0.4 27.8 28 

0.4 and 0.6 . . . . , 1 8.6 20 

0.6 and 0.8 12.2 12 

0.8 and 1 .0 $,6 4 

0.0 and 0.29 .... 49.8 50 
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The agreement shown here between the actual and theoretical 
percentages is somewhat closer than in the preceding case, the 
maximum discrepancy being only 1.6 and the average only 0.6 
per cent. 

Finally, the following data derived from one thousand act- 
ual errors may be cited. The errors of one hundred inter- 
polated values rounded to the nearest unit of the last tabular 
place were computed * for each of the interpolating factors 
0.1, 0.2, . . . 0.9. The averages of these several groups of act- 
ual errors are given along with the corresponding theoretical 
errors in the parallel columns below: 

Interpolating Actual Theoretical 

Factor. Average Error. Average Error. 

0.1 0.338 0.320 

0.2 0.288 0.303 

0.3 0.321 0.304 

0.4 0.268 0.290 

0.5 0.324 0.333 

0.6. . . 0.276 0.290 

0.7 0.321 0.304 

0.8 0.289 0.303 

0.9 0.347 0.320 

The average discrepancy between the actual and theoret- 
ical values shown here is 0.017. It is, perhaps, somewhat 
smaller than should be expected, since the computation of the 
actual errors to tliree places of decimals is hardly warranted 
by the assumption of dependence on first differences only. 

The average of the whole number of actual errors in this 
case is 0.308, which agrees to the same number of decimals 
with the average of the theoretical errors, f 

* By Prof. H. A. Howe. See Annals of Mathematics, Vol. Ill, p. 74. 
The theoretical averages were furnished to Prof. Howe by the author. 

f The reader who is acquainted with the elements of the method of least 
squares will find it instructive to apply that method to equation (i), Art. 13, 
and derive the probable error of e. This is frequently done without reserve by 
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Prob. 25. Apply formulas (3) of Art. 12 to the case of the sum 
or difference of two tabular logarithms and derive the correspond- 
ing values of the probable, mean, and average errors. The graph 
of 0(e) is in this case an isosceles triangle whose base, or axis of e, 
is 2, and whose altitude, or axis of 0(e), is i. 

those familiar with least squares. Thus, the probable error of ei or et being 
0.25, the probable error of e is found to be 



0.25 Vi — 2/-j- 2/*. 

This varies between 0.25 for /= o and 0.18 for ^ = ^ ; while the true value of 
the probable error, as shown by equations (3), Art. 13, varies from 0.25 to 0.15 
for the same values of /. It is, indeed, remarkable that the method of least 
squares, which admits infinite values for the actual errors ei and e^, should give 
so close an approximate formula as the above for the probable error of ^. 

Similarly, one accustomed to the method of least squares would be inclined 
to apply it to equation (4). Art. 13, to determine the probable error of e. The 
natural blunder in this case is to consider ^1, ^a, and <ra independent, and ^s like 
ei and e^ continuous betweer the limits 0.0 and 0.5 ; and to assign a probable 
error of 0.25 to each. In tMs manner the value 



0.25^2(1 -/ + /*) 

is derived. But this is absurd, since it gives 0.25 V2 instead of 0.25 for / = o. 
The formula fails then to give even approximate results except for values of / 
near 0.5. 
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